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1 Introduction and main results 

In this paper we construct discrete analogues of the Dixmier operators, that 
is, commuting difference operators corresponding to a spectral curve of genus 
1 whose coefficients are polynomials of the discrete variable n. 

The theories of commuting ordinary differential operators and commuting 
difference operators have a lot in common, therefore, to begin with, we recall 
some fundamental results concerning commuting differential operators. 

According to the Burchnall-Chaundy Lemma [1], if the two differential 
operators 

d n d n ^ d m d m ^ 

Ll = d^ +Un ~ l{x) d^ + - ' 12 = d^ +Vm - l{ ~ x) d^ + - ' - +V ° {X) 

commute, then there exists a nonzero polynomial Q(X, n) of two commuting 
variables A,/i, such that 

Q(L 1 ,L 2 ) = 0. 

On the plane C 2 with the coordinates (A, /i) this polynomial defines the spec- 
tral curve 

r = {(A,/i)eC 2 :g(A, /U ) = 0}. 

Suppose that T be a smooth curve. Considering some common eigen-function 
ip of the operators L x and L 2 , one can show that the corresponding eigen- 
values A and /i are connected by a polynomial relation Q(A,/i) = (see 
[2]). Thus, the spectral curve parametrizes the common eigen-functions and 
eigenvalues of L\ and L 2 

L 1 ij;(x,P) = \il>(x,P), L 2 il>(x,P)= t nj;(x,P), P=(X,fi)eT. 

The rank of operators Li, L 2 is called dimension of space of common eigen- 
functions with fixed eigenvalues (A, ji) G V. 



1 



In the case of operators of rank 1=1, one can find the function ip(x, P) 
(Baker-Akhiezer function) in explicit form by means of the theta-function of 
the Jacobi variety of the curve T, and in this case the operators' coefficients 
are easily found. Krichever [2] gives a classification of operators of rank / > 1. 
But the problem of finding the operators' coefficients is not solved in the 
general case. The following particular results relating to operators of rank 
more than 1 are known: Dixmier [3] found examples for operators of rank 2 
corresponding to an elliptic curve with polynomial coefficients 

where a is some constant. Krichever and Novikov [4] found all operators of 
rank 2 corresponding to the elliptic curve, and Grinevich [5] found spectral 
data corresponding to operators with rational coefficients. Mokhov [6] found 
operators for I = 3, also corresponding to the elliptic curve. In [7] and [8], 
we found examples for operators with rank 2 corresponding to the curve of 
genus (7 = 2, amongst which there are operators with polynomial coefficients, 
and also formally self-adjoined operators. 

As for smooth operators, for commuting difference operators of the form 

N+ M+ 

where n G Z is a discrete variable, and T — the shift operator on the discrete 
variable 

T/(n) = /(n + l), 

there exists a spectral curve T given in C 2 by a polynomial Q(X, //), parametri- 
zing their common eigen-functions and eigenvalues 

L 1 il>(n,P) = hl>{n,P), L 2 ^(n,P)=^(n,P), P=(\,fi)eT. 

Define the rank of operators Li and L 2 in the same way as for the smooth case 
as a dimension of space of common eigen-functions in the point in general 
position P G T. One of the basic differences between the discrete and the 
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smooth case consists of the following. Any commutative ring of ordinary 
differential operators is isomorphic to a ring of meromorphic functions on the 
algebraic curve with pole only in the point Q G T, whereas any commutative 
ring of difference operators is isomorphic to a ring of meromorphic functions 
on the algebraic curve with m poles, where m can be any natural number 
[9]. Such operators are called m-point operators. 

Mumford [10] and Krichever [11] found spectral data corresponding to 
two-pointed operators of rank 1. 

In this paper we will consider only single-pointed operators. For I > 1, 
finding the function ip(n, P) reduces to solving a Riemann problem, and it 
can not be found in the explicit form. Krichever and Novikov [9] (see also 
[4]) suggest a method for finding the operators' coefficients called method 
of Turin parameters deformation, and this method does not require to find 
ip(n, P). They show that the operators' coefficients can be restored from solu- 
tions of equations on Turin parameters of holomorphic stable vector bundles 
that are uniquely given by function ip(n,P), with the operators' coefficients 
depending on arbitrary / of functional parameters. Namely, they show that 
in order to restore the operators' coefficients it suffices to find the matrix 
function 

X (n, P) = *(n + 1, P)* _1 (n, P), 

where ^(n, P) is a Wronski matrix (see below) built on some basis in the 
space of common eigen-functions. Applying this method, Krichever and Novi- 
kov found operators of rank 2 corresponding to the elliptic curve. Simultane- 
ously, the operators' coefficients are expressed by ( and p-Weierstrass functi- 
ons of two functional parameters. 

In this paper we show spectral data for operators of rank 2 correspond- 
ing to the elliptic curve and whose coefficients are expressed by elementary 
functions of functional parameters. In particular, we find operators that have 
polynomial coefficients, just like the Dixmier operators. 

We can assume the affine part of the spectral curve V to be given in C 2 
with the coordinates (z, w) by equation 

w 2 = F(z) = + c 2 z 2 + dz + 1. 

Take Q — (0,1) e T as preferred point. Q — (0,1) e T (the ring of 
commuting difference operators will be isomorphic to a ring of meromorphic 
functions on T with pole in Q). The curve T allows a holomorphic involution 

a : T — > T, a(z, w) = (z, —w). 
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At / = 2 the matrice x( n > P)i P = (z,w) E T has the form 

X(n,P)= ( Xi (n,P) X2(n,P) 
We consider the case where the involution a does not change Xi, e.g. 

Xi(n,P)= X i(n,a(P)), PeT. (1) 
Then, the following holds: 

Theorem 1 The functions xi( x ,P) u X2(x,P) have the form 

Xl( „,j>) = _£<=L + c <"> 



z — 7(71) 7(71)— 7(71 + 1)' 
. 1 a(n) W7(n) ,, . 

p) = s + w^im + 2,^) - 7) + "<•*>• 

where 

7 (n-l)(a 2 H-F( 7 (n)) 
W 47(n)( 7 (n)- 7 (n-l)) ' 1 J 

(a(n + 1) - 1)7(71) + (a(w) + 1)7(71 + 1) 
W 2(7(71)- 7(71 + l)) 7 (n + l) 

7(71), a(n) are arbitrary functions of the discrete variable n G Z. 

The function Ai(z) on the curve T, with the only second-order pole in Q, 
looks as follows: 

1 Ci w 
1 = 2^2 + Tz + 2^2' 
Let bi(n) , ei(n) , Pi denote the series decomposition coefficients of functions 
Xi, X2, and Ai in the neighbourhood of Q: 

Xi(n, z) = 6o(n) + &i(n)z + . . . , X2(n, z) = - + e (n) + e x {n)z + . . . , 

z 

1 Pi c ? C 2 Ci 

The coefficients 6j,ej, i = 0, 1 are expressed by 7(71), and a{n) by formulae 
(6)-(9). 
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Corollary 1 The operator L(\i) has the following form 

L(Ai) = T 2 + Ul (n)T + uo(n) + Jt-iWr 1 + w_ 2 (n)T- 2 , 

Ml (n) = pi - e (n) - e (n + 1), 
tt (n) = p - b (n) - b (n + 1) - pie (n) + e„(n) - ei(n) - ei(n + 1), 

u_i(n) = -6i(n) + 6o(n) ^-pi - _ j| + e (n - 1) + e (n)^ , 

w_ 2 (n) = b {n)b (n - 1). 
The basic result of this work consists of the following: 
Theorem 2 // we put the function parameters a(n), 7(n) equal to 

a(n) =n + l, 7(71) = 7i, 
i/ien £/ie operators will have coefficients polynomial on n. At the same time, 

( n 4 1 1 

L 2 = T 2 + 2(n + 2)T - ( — + n 3 - -(1 - c 2 )n 2 - -(8 - a - c 2 )n 

~l(n 3 + (c 2 - l)n + ci - 2)(n 2 + n - 1)^+ 
( n3 + (C2 - l)n + ci - 2) (n 3 - 3n 2 + (2 + c 2 )n + Ci - c 2 - 2) (n + 1) (n - 2)T~ 2 . 

lb 

In Theorem 2 we do not list the third order operator due to its bulkiness. 
Let's look at an example. Let a spectral curve be given by the equation 



w 2 = z 4 + z 2 + 1, 



then 



L 2 = T 2 + 2(n + 2)T - - (n 4 + 2n 3 - 7n - 5) - - (n 3 - 2) (n 2 + n - l) T _1 

+ ^ (n 3 - 3n 2 + 3n - 3) (n + 1) (n - 2) T~ 2 , 
L 3 = T 3 + ^3n + y ^ T 2 - ^ (n 4 + 4n 3 + 5n 2 - 8ra - 14) T 



~ (2n 5 + 7n 4 + 10n 3 + n 2 - Yin - 5) 

+ — (n 8 - 2n 6 - 12n 5 - 3n 4 + 10n 3 + 10n 2 + 6n - 12) T _1 
16 

+— n {in 2 - n - 5) (n 6 - 3n 5 + 3n 4 - 5n 3 + 6n 2 - 6n + 6) T~ 2 — 

— (n - 3) (n 2 - l) (n 3 - 2) (n 3 - 6n 2 + 12n - 10) (n 3 - 3n 2 + 3n - 3) T" 3 . 

In [7] and [8], we use an analogue of reduction (1) in the smooth case in 
order to find commuting differential operators of rank 2, genus 2. 

In section 2 we recall the Krichever-Novikov equation [9] on the discrete 
Turin parameters dynamics and give formulae for coefficients of operators of 
rank 2, genus 1 that are expressed by ( and jp-Weierstrass functions. 

In section 3, theorems 1 and 2 are proven. 

2 Krichever-Novikov equations on the discrete 
dynamic of the Turin parameters 

As noted above, for I > 1 function if)(n, P) cannot be found in its explicit 
form. Let ty(n,P) denote a Wronski matrice ^(n,P) = ip j (n + i), where 
^(n, P), 1 < j < I is some basis in the space of common eigen-functions. As 
shown in [9], the number of zeros of functions det^(n, P) equals Ig, where 
g is the genus of T. Denote them by 7i(n), . . . , 7j 9 (n). By aij(n) we denote 
vectors where 

These vectors are defined up to proportionality. The set (jjin), (Xj(n)) is 
called Turin parameters. The Turin parameters uniquely define a stable holo- 
morphic bundle of rank / on the curve T. The following theorem, proven in 
[9], defines the discrete dynamics of Turin parameters. 

• The matrix function x(n, P) has simple poles in points 7j(?i). The re- 
lations on the residues of the matrix elements 

^Res 7s(n)X ™(n, P) = a\{n)Res ls{n)X m ' f (n, P) (3) 

do hold. 



6 



The points 7s (n + 1) are zeros of the matrix determinant P), e -g- 

det X (n, 7s (n + l)) = 0. (4) 
The vector aij(n + 1) satisfies the equation 

a J -(n + l)x(n, 7j (n + l)) = 0. (5) 

Except for the poles jj(n) one of the components x has a simple pole in 
the preferred point Q. 

As follows from the Riemann-Roch theorem, x( n -> P) is uniquely restored 
by the Turin parameters ( 7 j(n), otj{n)) and by / arbitrary functional parame- 
ters. 

This theorem enables us to find operators of rank 2 corresponding to an 
elliptic curve. Below, we give the according theorem, for simplicity only for 
the case where a holomorphic involution on an elliptic curve switches the 
location of the poles of function Xi- We give the elliptic curve T as a factor 
C/A, where A is some lattice in C and let z be a coordinate in C. Assume 
that points 71 (n) and 72 (n) are interchanged by a holomorphic involution on 

r 

°{li{n)) = 72 W, 
where a(z) = —z. The following theorem was proven in [9]. 

• The operator corresponding to function p(z) has the form 

L = L\- p( 7 (n)) - p( 7 (n - 1)), 

where 

L 2 =T + v(n) + c(n)T-\ 
c(n + 1) = - A (s 2 (n) - l)F( 7 (n + 1), 7 (n))F( 7 (n - 1), 7 (n)), 

v(n + 1) = i(s(n)F( 7 (n + 1), 7 (n)) - s(n + l)F( 7 (n), 7 (n + 1)), 

u) = C(« + «) - C(« - - 2C(u), 
s(n), 7 (n) are arbitrary functional parameters. 

As can be seen from the formulae for the coefficients of operator L, the 
coefficients are expressed by function parameters s(n) and 7 (n) applying 
the p and (- Weierstrass functions. In corollary 1 and theorem 2 we found 
conditions for these coefficients being elementary functions and, in particular, 
polynomials. 
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3 Proof of theorems 1 and 2 

As follows from our assumption on the invariance Xi{ n -> P)i P — (z,w) & F 
under the action of involution a, Xi( n , P) has the form 

c(n) , , . 
Xi(n,P)= y +do(n), 
z — 7 (n) 

where c(n), 7 (n) and do(n) are some functions of a discrete variable. Note 
that function xi( n ,P) has poles in 

P{n) = ( 7 (n), VP(l(n)), aP{n) = ( 7 (n), -y/F^n)). 

From equalities (4), being equivalent in our case to 

X i(n, Pin + 1)) = X i(n, <xP(n + 1)) = 

get 

. c(n) dn) 



z — 7 (n) 7 (n) — 7 (n + 1) 
Look for function X2 in the form 

. 1 a(n) wy(n) ,, . 

P) = T; + 2(7^)) + 2»( 7 («) - 7) + d(n) ' 

where a(n) and d(n) are some functions. Note that x 2 (n, P) has poles in 
points Pin),aPin) and Q, moreover ResgX2 = 1- 

Vectors «i and 0:2 are defined up to proportionality, therefore we can 
assume them to have the form 

ai(n) = (ai(n),l), a 2 (n) = (a 2 (n), 1). 

Then from (5) we get 



2 7 (n) 2( 7 (n) - 7 (n - 1)) 2 7 (n)( 7 (n - 1) - 7 (n)) 
a 2 (n) = -X 2 (n - 1, oP(n)) 



8 



1 + , „- + »-l). 



2 7 (n) 2(7(11) - 7(11 - 1)) 2 7 (n)(7(n - 1) - 7(11)) 

Further, the residues of functions Xi an d X2 in poles P(n) and aP(n) are 
equal to 

Res P(n )Xi = Res CTP(n )Xi = c(n), 
Res P(n) X2 = ^(a(n) - ^(7(71)), 

Res a p (n) X2 = ^(a(n) + y/F^n)). 
From equalities (3), in our case taking the form 

Res P(n) Xi = a(n)Res P(n) X2, 
Res CT p (n) Xi = a(n)Res CT p (n) X2, 



get 



d(n) 



7 (n-l)(a 2 (n)-P( 7 (n)) 
° {n) 4 7 (n)( 7 (n)- 7 (n-l) ' 

(a(n + 1) - 1)7(71) + (a(n) + l) 7 (n + 1) 



2( 7 (n) - 7 (n + l)) 7 (n + l) 

Thus, theorem 1 is proven. 

Now find the coefficients of operator L 2 as follows. Express ip{n + 2, P) 
and ip{n — 2, P) by ip{n — 1, P), ^(n, P), xi and X2- For this use the identity 

ip(n + 1) = ip(n - l)xi(n) + ip(n)x2(n) 

from which 

ip(n + 2) = tp(n - l)xi(n)x2(n + 1) + ^(n)(xi(n + 1) + * 2 Wx 2 (w + 1)), 

ip(n - 2 = -^(n - 1 — + — -. 

Xi(n-l) xi(^-l) 

Now replace T 2 ip(n) and T~ 2 ip{n) in equation 

L 2 ^ = (T 2 + Ul (n)T + tio(n) + M-iT" 1 + u_ 2 T- 2 )^(n) = A^(n) 
by corresponding expressions. Get 

Pi(n, P)^(n, P) + P 2 (n, P)^(n - 1, P) = Ai^(n, P), 
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where 

u- 2 (n) 



Pi( n ) = Xi( n + 1) + Xi{ n + l)X2(n) + u 1 (n)x2(n) + u (n) + 
Pz(n) = X2(n + l)xi(ri) + ui(n)xi(n) + u_i(ra) - w_ 2 (n 



Xi(n- 1)' 
X2(n- 1) 



Xi(n- 1)' 

In the space of common eigen-functions of operators L 2 and L 3 we can chose 
a basis and t/> 2 normalized by conditions 

Mno,P) = 1, ^ 2 (n ,P) = 0, 

with functions Xi)X2 not depending on normalization points no (see [9]), 
consequently, we have the identities 

P 1 (n,P)-A 1 = 0, P 2 (n,P) = 0. 

Now, having found in theorem 1 the functions Xi and X2, find the coefficients 
Ui(n) of operator L 2 from the series decomposition P\ — Ai and P 2 in the 
neighbourhood of point Q. 

The series decomposition coefficients bj,ej of functions Xi an d X2 in the 
neighbourhood of Q are expressed by the functional parameters a{n) and 
7(71) according to formulae 

7 (n-l) 7 (n+l)(F( 7 (n))-a 2 H) 
01 J 4( 7 ( n -l)-7( n ))(7( n )-7(n + l)) 7 («) 2 ' {) 

= 7 (n-l)(a 2 (n)-F( 7 (n)) 
U J 4( 7 (n- 1) - 7 (n)) 7 3 (n) ' 1 1 

. 1 ( c\ 1 a(n) 
2 V 2 7 (n) 7(71) 



(o(n + 1) - l)7(n) + (o(n) + 1)7(71 + 1] 
(7(71) - 7(71 + l))7(n + 1) 



(8) 



8 - 8a(w) + 4cg(n) - (c\ - 4c 2 ) 7 2 (n) 

ei(n) = le^H ' (9) 

Analogically, find the operator corresponding to a meromorphic function with 
one third order pole in Q. Direct verification yields that for 

a(n) =n + l, 7(71) = n 
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the operators' coefficients are polynomials on n. Thus, theorem 2 is proven. 
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